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Abstract. By analyzing the Rossi X-ray Timing Explorer data on SAX J1808.4— 3658, we show that the X-ray emission 
in accretion powered millisecond pulsars can be produced by Comptonization in a hot slab (radiative shock) of Thomson 
optical depth Tes ~ 1 at the neutron star surface. The escaping radiation consists of two components: a black body and 
a hard Comptonized tail. These components have very different angular distribution: the black body peaks along the slab 
normal (a "pencil"-like emission pattern), while the tail has a broader angular distribution (a "fan"-like pattern). This results 
in very different variability properties. We construct a detailed model of the X-ray production accounting for the Doppler 
boosting, relativistic aberration and gravitational light bending. We are able to reproduce the pulse profiles at different energies, 
corresponding phase lags, as well as the time-averaged spectrum. We obtain constraints on the neutron star radius: i \ km 
if its mass M = 1.6M,;,, and ~ 8.5 km if M = 1.4M0. We present simple analytical formulae for computing the light curves 
and oscillation amplitudes expected from hot spots in X-ray bursters and accretion powered millisecond pulsars. We also 
propose an analytical expression that can be used to determine the size of the black body emission region from the observed 
properties. 



INTRODUCTION 

In order to understanding the X-ray production in accret- 
ing black hole and neutron star sources in X-ray binaries, 
one needs to model their spectral and temporal prop- 
erties simultaneously. In the black hole case, the emis- 
sion is believed to originate in the accretion disk and/or 
its corona. Detailed spectral fits to the broad band spec- 
tra leave little doubts that the main emission mechanism 
is thermal or non-thermal Comptonization. However, in 
spite of extensive efforts, the exact geometry of the emis- 
sion region is still not known 1 12]. The hard X-rays can 
be produced in the central hot flow flTil and/or in the 
magnetically dominating corona/jet in the vicinity of a 
cooler accretion disk 111 llOll . Temporal variability pro- 
vide some additional constraints on the size of the emis- 
sion region |9], but has difficulties in breaking the model 
degeneracy 1 13]. The observed spectral-temporal corre- 
lations seem, however, to favor the hot inner disk sce- 
nario 0, 0, |23l at least for the hard state dominated by 
thermal Comptonization. 

In the case of accreting neutron stars in low-mass X- 
ray binaries, the emission can be produced in the ac- 
cretion disk (and corona) as well as at the neutron star 
surface, with the later dominating the total energy out- 
put ]19, 22]. Depending on the magnetic field strength 
and the accretion rate, there can be different scenarios. 
The gravitational energy can be dissipated gradually in 
the boundary (or spreading) layer, while the matter is 



slowing down from the rapid Keplerian rotation to a 
slower rotation of a neutron star, or abruptly in the shock 
when the accreting material, following magnetic field 
lines, is channeled towards magnetic poles of the star. Di- 
rect spectral decomposition can be quite complicated jist]. 
The analysis of the Fourier frequency resolved spectra of 
neutron stars in low-mass X-ray binaries shows that the 
accretion disk producing soft X-ray emission does not 
vary much, while the hard radiation (from the boundary 
layer?) is responsible for most of the observed variability 
at high Fourier frequencies ] 8 ] . 

In all the cases considered above, it is not absolutely 
clear what is the origin of a given spectral component. 
In accretion powered ]24] as well as nuclear powered 
(X-ray bursters) ]20] millisecond pulsars, the situation 
seems to be simpler Observations of the coherent pulsa- 
tions confirms that the bulk of the emission is produced at 
the surface of a neutron star which rotates at 200-600 Hz. 
We can also fold the light curve on a pulsar period and 
study the pulse profile as a function of energy. For accret- 
ing pulsars the statistics is better since one can fold the 
profiles over a longer observational period (days rather 
than seconds in the X-ray burst oscillations). Detailed 
modeling of the pulse profiles allows us to get interesting 
constraints on the compactness of the neutron star as well 
as on the emission pattern from the neutron star surface 
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FIGURE 1. Time-averaged spectrum of SAX J1808.4-3658 
as observed by RXTE and the model spectrum. See text for 
details. 

SPECTRUM AND RADIATION PATTERN 

We assume that the emission originates close to the neu- 
tron star surface in an accretion shock which we ap- 
proximate as a plane-parallel slab. The time-averaged 
observed spectrum of the best studied accretion pow- 
ered pulsar SAX J1808.4-3658 (see Fig. [Hi is similar 
to those in other neutron star and black hole sources. It 
can be represented as a composition of a black body with 
temperature fcTbb = 0.66 keV (dash-dotted curve), ther- 
mal Comptonization (dots) by electrons of temperature 
^ 50 — 100 keV and optical depth Tgs ^ 1, Compton re- 
flection bump and the iron line (dashes). The Compton 
reflection (from the disk) has a rather small amplitude 
because of the disruption of the accretion flow in the neu- 
tron star vicinity f3\ . 

If the settling bulk velocity smaller than the thermal 
electron velocity, we can easily compute the angular dis- 
tribution of the escaping radiation jl6ll23ll . The radiation 
pattern is very different for photons of different scattered 
orders (see Fig.|2}. One sees that the black body (marked 
with 0) is strongly peaked along the normal to the slab, 
while photons scattered many times are beamed in the di- 
rection making an angle ^ 50° — 60° from the normal. In 
the Comptonization process, scatterings also shift pho- 
tons along the energy axis so that higher energy photons 
are scattered more times. 



MODEL AND RESULTS 

The light curves from a circular spot are computed 
accounting for the special relativistic effects (Doppler 
boosting, relativistic aberration) and the gravitational 
light bending in Schwarzschild metric. The model is 
described in details in lITsll . The spot spectrum is as- 




FIGURE 2. Polar diagram of the normalized radiation flux 
escaping from an electron scattering slab of Thomson 
optical depth Tes — 0.7. Here arccos fl is the angle between the 
normal and the photon direction. Different scattering orders are 
shown and marked by numbers. 
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FIGURE 3. Light curves expected from a slowly rotating 
star (dashed curves) and that rotating at 401 Hz (solid curves). 
Parameters: M = IAMq, R = 2rg = 8.4 km, i = 80°, = 11°, 
T = 0.16 and a = —0.78 (this parameter corresponds to the 
scattering optical depth Tgs ~ 0.7). From IT^ . 
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FIGURE 4. (a) Pulse profiles of SAX J1808.4-3658 in 
the 3-4 keV (circles) and 12-18 keV (squares) energy band 
and the model light curves (solid and dashed curves). Same 
parameters as in Fig. |3| (b) The angular distribution of the 
intrinsic black body (solid curve) and Comptonized (dashed 
curve) fluxes /x/(/i) in the spot co-rotating frame normalized 
as / fll{il)dfl — 1. Only the range of angles between the dotted 
lines is actually observed, (c) The observed (crosses) and the 
model (solid curve) phase lags at the pulsar frequency relative 
to the 3-4 keV band. From 115,1 . 



sumed to consist of two components: a black body and a 
Comptonization tail with different angular distributions. 
Their spectral shapes are taken as in the observed time- 
averaged spectrum (see Fig. 0. The model parameters 
are: the pulsar frequency (fixed at v = 401 Hz for SAX 
J1808.4— 3658), the neutron star mass M, stellar radius 
R, inclination /, colatitude of the spot center 9, "optical 
depth" T, that describes the angular dependence of the 
black body intensity /bb(M) °= 6xp(— T//i), and a param- 
eter a determining the angular distribution of the Comp- 
tonized radiation /sc(/i) °^ 1 +fl/X. This linear dependence 
of the intensity mimics the angular distribution of radia- 
tion escaping from a slab (Fig.|5}- 

A rapid rotation of the star causes significant changes 
in the pulse profile. We illustrate this by showing in 
Fig.|3lthe light curves of a slowly rotating star (dashed 
curves) and those modified by the Doppler boosting and 
aberration for a neutron star of rotational frequency 401 
Hz (solid curves). We consider both the black body and 
the Comptonized emission. The pulse profiles strongly 
depend on the assumed emission pattern. Light bending 
depends on the compactness of the star R/rg (where 
rg = 2GM/c^) and reduces the variability amplitude with 
respect to a less compact star (see below). 
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FIGURE 5. Constraints on the radius of the compact star 
from the light curve of SAX J1808.4-3658 are shown by 
circles with double error bars (corresponding to 90% and 99% 
confidence limits). Different equations of state for strange and 
neutron stars are shown for comparison. See IT^ for details. 



The model parameters can be constrained by fitting 
the observed light curves from SAX J1808.4-3658 (see 
Fig. 13. The fitted emission pattern of Comptonized and 
black body radiation confirms our expectations: the black 
body is beamed along the surface normal while the hard 
radiation flux (proportional to the intensity times the 
cosine of the projected area) is more isotropic (Fig.|^). 
One of the main result of this study is determination 
of the radius of the compact star The dependence of 
the radius on the assumed mass is shown in Fig.|5] For 
a standard neutron star mass of 1.4M0, we get R = 
8 .5 ± 0.5 km, while for a more massive star M = 1 .6M0, 
the radius of 11 km is consistent with some neutron star 
equation of state. The inclination of the system can be 
constrained to / > 65°. Our model is able to fit the time- 
averaged spectra, the energy dependent pulse profiles and 
the observed phase lags, simultaneously. 

Very similar pulse profiles (but with higher oscillation 
amplitude) were also observed from the recently discov- 
ered fifth ms pulsar XTE J1814-338 |21]. In the frame- 
work of the present model, the data can be explained for 
example, by increasing the colatitude of the spot center 9 
from 11° to ^ 17°. The presence of the harmonic in the 
burst oscillations \2]\ could result from the larger contri- 
bution of the Comptonized radiation (with a "fan"-like 
emission pattern) to the observed flux in this source com- 
paring to other X-ray bursters since a black body does 
not produce harmonics (see eq.|5}- This interpretation is 
supported by the observed hard spectrum flUl . 



ANALYTICAL LIGHT CURVES, 
OSCILLATION AMPLITUDES, AND 
SPOT SIZE 

Using formalism described in |01 (see fisll for details) 
we obtain simple expressions for the light curve from a 
small spot at a rapidly rotating neutron star. For example, 
the bolometric flux from a black body spot is 



dF=(l--Al 8 
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— + 1 - — cos 1/ 



R 



dS 

^0^, (1) 



where /q is the black body intensity (in the spot comoving 
frame), dS is the spot area, D is the distance to the 
source, cosy/ — cosi cos + sin/ sin0 cos0, is the 
pulse phase, and 5 is the Doppler factor. We can also get 
simple (but very accurate) expression for the oscillation 
amplitude (here from one spot and neglecting Doppler 
boosting which does not affect it much): A = (/vnax — 
^min)/(fmax +-Pmin) = U /Q, where we defined t/ = (1 - 
rg/R) sin/ sin and Q = / R + {I ~ rg/R) cos i cos 9 . 

The time-av erag ed flux observed from a spot of angu- 
lar radius p is llSll (assuming the spot is always visible) 



2 r 



e+^tan2^ 
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kTtR^sin^p/D^. (2) 



This formula allows us to obtain a simple expression for 
estimating the emission region size (compare to 1 18]): 



R sin p = roc 



e + ^tan^^ 
^ R 2 



-1/2 



(3) 



which can be solved by iterations. Here Too is the "ob- 
served" spot size, = <y&BTti~L/D^^ and Too is the fitted 
temperature. For small p, we get Rp = TooQ^^I^. For ex- 
ample, if we look at the spot always along the symmetry 
axis (take / = 0, Q — 0), we get Q~ \ and Rp = r„o. Thus 
in that case, light bending does not affect the size. A spot 
of a finite size produces a smaller variability amplitude 
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For an anisotropic source given by a linear relation 
/(jU) o<: 1 + fl/X at a slowly rotating star (5 — 1), we 
can obtain the ratio of the variability amplitudes at the 
fundamental frequency Aq and the harmonic A\, using 
equation ([l): 

Ai _ aU/2 



(5) 



1 +2aQ 

One sees that a black body source a = does not produce 
any harmonics. We can estimate the impact of rapid rota- 
tion on that ratio. Expanding Doppler factor in powers of 
the equatorial velocity /3eq = Veq/c and keeping the first 
term only, we get (see Iil4il for details): 



j3eq\/l — rg/R sin/ sin0. 



(6) 



Thus this ratio depends linearly on sin/ and sin0 (as 
for anisotropic source, eq. |5}. One also sees that due 
to a linear dependence on the rotational frequency, it 
is much easier to get large amplitude of the harmonic 
with an anisotropic source than just by rotation. These 
expressions can be used for analysis of the data from ms 
pulsars and X-ray bursters. 
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